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Abstract 

The Asteraceae family is at the forefront of the evolution due to frequent hybridization. Hybridization is associated with the induction 
of widespread genetic and epigenetic changes and has played an important role in the evolution of many plant taxa. We attempted 
the intergeneric cross Chrysanthemum morifolium x Leucanthemum paludosum. To obtain the success in cross, we have to turn to 
ovule rescue. DNA profiling of the amphihaploid and amphidiploid was investigated using amplified fragment length polymorphism, 
sequence-related amplified polymorphism, start codon targeted polymorphism, and methylation-sensitive amplification polymor- 
phism (MSAP). Hybridization induced rapid changes at the genetic and the epigenetic levels. The genetic changes mainly involved loss 
of parental fragments and gaining of novel fragments, and some eliminated sequences possibly from the noncoding region of 
L paludosum. The MSAP analysis indicated that the level of DNA methylation was lower in the amphiploid (-45%) than in the 
parental lines (51 .5-50.6%), whereas it increased after amphidiploid formation. Events associated with intergeneric genomic shock 
were a feature of C. morifolium x L paludosum hybrid, given that the genetic relationship between the parental species is relatively 
distant. Our results provide genetic and epigenetic evidence for understanding genomic shock in wide crosses between species in 
Asteraceae and suggest a need to expand our current evolutionary framework to encompass a genetic/epigenetic dimension when 
seeking to understand wide crosses. 
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Introduction 

The whole-genome doubling of natural wide hybrids to pro- 
duce allopolyploids is a common event in the plant kingdom, 
and their formation has provided a significant means of 
promoting adaptation and speciation (Kearney 2005; Rapp 
and Wendel 2005). Plant genomes also have the capacity to 
spontaneously double their genomes, which in turn produce 
autopolyploids, whereas autopolyploids are less frequently en- 
countered than allopolyploids (Grant-Downton and Dickinson 
2006; Paun et al. 2006). Current estimates suggest that at 
least 70% of modern plant species are either polyploids or de- 
scendants from a polyploid ancestor. Polyploid genomes are 
fundamentally unstable, tending to decay back to the 



functional diploid state through a process of pseudogeniza- 
tion and segmental deletion known as diploidization (Wolfe 
2001; Ma and Gustafson 2005). Arabidopsis, maize, and rice, 
traditionally regarded as diploid species, are all good examples 
of diploidized polyploids (Chen 2007). The first few genera- 
tions following a wide hybridization event, as has been dem- 
onstrated for a diversity of synthetic allopolyploids, are prone 
to genomic perturbation at both genetic and epigenetic levels, 
a phenomenon referred to as "genomic shock" (Yu et al. 
2005; Jackson and Chen 2010). 

The nature of genomic shock has been repeatedly investi- 
gated (Matsuoka 2011). Its manifestation includes 
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chromosomal rearrangement, the gain and loss of chromo- 
some segments, gene repression and activation, subfunctio- 
nalization, and transposon activation, as well as changes in the 
epigenome, in particular with respect to patterns of cytosine 
methylation (Adams et al. 2003; Paun et al. 2006; Otto 2007; 
Kawakami et al. 2011). Nevertheless, the plant seems well 
tolerated in many genomic groups by combining traits from 
two parental species and even potentially benefits from 
those changes, which promotes growth vigor in response to 
"genomic shock" stress in their hybrids (McClintock 1984; 
Paun et al. 2007; Ungerer and Kawakami 2013). Some of 
these alterations occur very rapidly after the hybridization 
and/or polyploidization event, as shown in wheat and 
Tragopogon (Shaked et al. 2001; Tate et al. 2006). Changes 
in global patterns of cytosine methylation, as detected 
by methylation-sensitive amplification polymorphism, have 
been directly associated with the hybridization and/or poly- 
ploidization event (Shaked et al. 2001; Salmon et al. 2005; 
Qi et al. 2010; Hegarty et al. 2011). However, there is little 
evidence for any rapid methylation changes in polyploidy for- 
mation of cotton and ploidy watermelon (Liu et al. 2001; 
Wang et al. 2009). Thus, polyploidy in the plant kingdom 
may involve an array of molecular events that are dependent 
on the genomic context. 

In general, the more distantly related are the prospective 
parents of a wide cross, the harder the cross is to make. In 
the context of crop improvement, wide crosses such as 
those involving parents from different genera or/and differ- 
ent ploidy levels are of particular interest because they 
offer the potentiality of introgressing novel genes into the 
primary gene pool of the crop (Ozkan et al. 2001; Zhang 
et al. 2011). The Asteraceae family comprises more than 
25,000 geographically widely dispersed species, represent- 
ing 10% of all angiosperms, and many of these species are 
polyploid (Bremer and Humphries 1993; Guo et al. 2012). 
Consequently, the Asteraceae family could provide a 
unique opportunity to understand the evolutionary geno- 
mics of lineage radiation and diversification at various phy- 
logenetic scales (Kane et al. 2011). Although transcriptome 
shock and cytosine methylation have been studied in inter- 
specific hybrids in Asteraceae (Hegarty et al. 2006; Tate 
et al. 2006; Hegarty et al. 2008), genetic and epigenetic 
alterations under intergeneric genomic shock are still poorly 
understood. 

The ornamental species chrysanthemum (Chrysanthemum 
morifolium) is a hexaploid, whereas the ploidy level of other 
Chrysanthemum spp. varies from diploid to decaploid (Kondo 
et al. 1 999). The perennial Asteraceae species Leucanthemum 
paludosum (formerly Mauranthemum paludosum), a native 
diploid species of the Mediterranean region, is also used for 
ornamental purposes. Here, we describe the producing of an 
allopolyploid from the cross C. morifolium x L. paludosum 
using ovule rescue and detail some of the major genetic and 
epigenetic changes induced by this intergeneric hybridization 



and polyploidization using a variety of DNA-based marker 
technologies. 

Materials and Methods 

Plant Material, Crossing, and Embryo Rescue 

The plant materials used are the C. morifolium x L. paludo- 
sum hybrid, chromosome-doubled hybrids derived from 
the cross combination, its parents C. morifolium cv. 
"Zhongshanzigui," and an accession of L paludosum (fig. 
^A-D) and the diploid species C. zawadskii. The materials 
are all maintained by the Chrysanthemum Germplasm 
Resource Preserving Centre, Nanjing Agricultural University, 
China (32°05'N, 118°8'E, 58 m altitude). The plants were 
propagated by cuttings. All plants were grown in a green- 
house held at 22 °C during the day and above 15°C during 
the night, with a relative humidity of 70-75% and under 
natural light. To make the wide cross, at 9:00-10:00 AM 
on a sunny day, the bisexual tubular florets of 
"Zhongshanzigui" were removed using tweezers before an- 
thesis and the inflorescences were enclosed within a paper 
bag when the stigmas first became visible on the developing 
receptacle. After 2 or 3 days, fresh pollen from newly blos- 
soming flowers of the L. paludosum was brushed onto the 
"Zhongshanzigui" pistil (at a Y-shaped stigma stage) using a 
fine-haired paintbrush. After pollination, the flowers were 
immediately bagged. Four, 6, 8, 10, 12, 15, 18, and 21 
days after pollination (DAP), individual pollinated florets 
were collected, fixed in FAA (1:1:18, formalin:glacial acetic 
acid:70% v/v ethanol), and stored at room temperature. 
Before examination, the florets were dehydrated by passing 
through an alcohol series, infiltrated with xylene, and em- 
bedded in paraffin wax, following the methods described by 
Deng et al. (2010). Sections of thickness 6-10jim were cut, 
stained in Heidenhain's hematoxylin, and then observed by 
light microscopy. For ovule rescue, ovaries were excised be- 
tween 10 and 15 DAP, surface sterilized by immersion in 
75% v/v ethanol for 35 s, followed by 10% v/v H 2 0 2 for 
15min, and rinsed six times in sterile water. The ovary coats 
were aseptically removed to extract the ovules; the extracted 
ovules were then placed on solidified MS medium 
(Murashige and Skoog 1962) containing 2mg/l 6-benzyla- 
denine (6-BA) and 0.5mg/l ot-naphthaleneacetic acid 
(Tang et al. 2009) and cultured at 25 °C under a 16-h 
photoperiod provided by cool-white fluorescent lamps 
(30jimol rrT 2 s~ 1 ). Regenerating plantlets were transferred 
to fresh solidified MS medium supplemented with 1.0% 
w/v sucrose (pH 5.7) until they had developed roots 1 cm 
long (Murashige and Skoog 1962), and the humidity was 
gradually lowered over a period of 1 month, after which 
they were potted into a 2:2:1 mixture of perlite, vermiculite, 
and leaf mold, and grown in a greenhouse. 
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Fig. 1. — Morphology of materials. (A) Chrysanthemum morifolium "Zhongshanzigui," (B) amphihaploid, (0 amphidoploid, and (D) Leucanthemum 
paludosum). From left to right of each line: plant morphology; floral morphology, bar: 1 cm; ligulate flower, bar: 2 mm; tubular flower, bar: 2 mm; leaf: 
0.5 cm. 



Cytological Confirmation and Genomic In Situ 
Hybridization Analysis 

Young root tips collected from parental cuttings and from the 
putative amphihaploid plantlets were held in ice water for 
20-24 h, then fixed in 3:1 ethanohglacial acetic acid (v/v) 
at 4°C for 24 h, and squashed in a drop of 45% acetic 
acid. Mitotic chromosome spreads were observed by phase 
contrast microscopy. The mitotic chromosomes of the puta- 
tive amphihaploid were also subjected to genomic in situ 



hybridization (GISH) by probing with genomic DNA of C. 
zawadskii extracted from young leaves using the CTAB 
(cetyltrimethylammonium bromide) method (Stewart and 
Via 1993) and labeled with biotinylated 16-dUTP (Roche 
Ltd.). The GISH procedure followed Deng et al. (2010). 

Chromosome Doubling 

Following the procedure described by Liu et al. (201 1), nodal 
segments from 1 -month-old hybrid plantlets were immersed 



Genome Biol. Evol. 6(1):247-259. doi:10.1093/gbe/evu008 Advance Access publication January 9, 2014 



249 



Wang et al. 



GBE 



in 600 mg/l colchicine for 48 h, then rinsed three times in ster- 
ile water, and placed on hormone-free MS medium for 
30 days, at which time the developed lateral buds were trans- 
ferred onto a rooting medium (Liu et al. 201 1). 

Genotypic Analysis 

DNA was extracted from the fourth and fifth leaves of three 
individual plants of the two parental lines, the three amphiha- 
ploids, and the corresponding amphidiploids using a Nuclei 
Isolation kit (Solarbio, China). The DNAs were subjected to 
amplified fragment length polymorphism (AFLP) profiling as 
described by Vos et al. (1995). Five hundred nanograms of 
DNA was digested overnight at 37 °C with 1 0 U each of EcoRI 
and Msel (New England Biolabs, Beijing, China). Then adap- 
tors were ligated to the ends of the digested fragments at 
16°C for 4h; the ligation mix contains the digested frag- 
ments, 5-pmol EcoRI and 50-pmol Msel adaptors (supplemen- 
tary table S1, Supplementary Material online), and 4 U T4 
DNA polymerase (NEB), after which the enzymes were inacti- 
vated by heating to 70 °C for 1 5 min. The amplicon generated 
from this template based on EcoRI and Msel primers lacking 
any selective bases (supplementary table S1, Supplementary 
Material online) was diluted 1:30 in ddH 2 0 and used as the 
template for subsequent selective amplifications based on 
EcoRI and Msel primers carrying three selective bases; in 
total, the nine primer combinations are EcoRI selective 
primer #2 plus Msel selective primer #5 (abbreviated 
"E2 + M5"), E2 + M6, E3 + M2, E4+M3, E4+M8, E6 + M7, 
E7 + M3, E8+M3, and E8+M7 (supplementary table S1, 
Supplementary Material online). The EcoRI primers were la- 
beled with 5-FAM (Carboxyfluorescein). Each polymerase 
chain reaction (PCR) and electrophoretic separation was re- 
peated twice. The amplicons were separated using an ABI 377 
DNA Analyzer device, and only fragments ranging from 80 to 
450 bp were considered. 

The same genomic DNAs were also subjected to sequence- 
related amplified polymorphism (SRAP) analysis, a DNA finger- 
printing technique that targets open reading frames (Li and 
Quiros 2001 ). In all, 24 SRAP primer combinations were used, 
involving 7 forward and 8 reverse primers (supplementary 
table S1, Supplementary Material online). The combinations 
wereME-1 combined with EM-1 (abbreviated "ME1 + EM1"), 
ME1+EM7, ME1+EM14, ME3 + EM2, ME3 + EM5, 
ME3 + EM10, ME6+EM1, ME6 + EM2, ME6+EM5, 
ME6+EM14, ME10+EM1, ME10 + EM7, ME10+EM11, 
ME10+EM14, ME14+EM2, ME14+EM5, ME16+EM2, 
ME16+EM10, ME16+EM11, ME16 + EM14, ME17 + EM1, 
ME17 + EM5, ME17 + EM7, and ME17 + EM15. Each 25-|il 
reaction mix comprised 20 ng genomic DNA, 200 jiM dNTP, 
1 .5 mM MgCI 2 , 0.6 jiM of each primer, 2.5 jil 1 0x buffer, and 
2 U Taq polymerase (Takara). The reactions were subjected to 
a 5-min denaturation at 94°C, then five cycles of 94 °C/1 min, 
35°C/1min, 72°C/2min, followed by 35 cycles of 



94°C/1 min, 50°C/1 min, 72°C/2min, and finally an exten- 
sion step of 72°C/7min. For the start codon-targeted poly- 
morphism (SCoT) analysis, the primers (supplementary table 
S1, Supplementary Material online) were derived from a con- 
sensus sequence in the flanking regions around the ATG start 
codon (Joshi et al. 1997). Each 25-jil PCR reaction mix con- 
tained 20 ng genomic DNA, 200 jiM dNTP, 1.5 mM MgCI 2/ 
0.6 jiM of each primer, 2.5 jil 10x buffer, and 2 U Taq poly- 
merase (Takara, Japan), and the cycling regime consisted of an 
initial denaturation step (94°C/3min), followed by 35 cycles 
of 94°C/1 min, 50°C/1 min, 72°C/2min, and a final exten- 
sion step of 72°C/5min. Both the SRAP and SCoT 
amplicons were electrophoresed through 8% nondenaturing 
polyacrylamide gels run at 300 V for 2.5 h in 1 x TBE buffer 
and were visualized by silver staining. 

Isolating and Sequencing Informative AFLP and 
SRAP Fragments 

Informative AFLP amplicons were also separated by electro- 
phoresis through 6% denaturing polyacrylamide gels and silver 
stained (Vos et al. 1995; Chen et al. 2007). Polymorphic AFLP 
fragments and polymorphic SRAP fragments were cut from 
the gel and eluted in 50 jil ddH 2 0 by boiling for 1 0 min. After 
centrifugation (1 2,000 x g, 5 min), an aliquot of the superna- 
tant was collected and served as the template. The fragment 
was reamplified using the same primer combination and PCR 
conditions. The resulting amplicon was inserted into the clon- 
ing vector PMD19 TA (Takara) and sequenced. The similarity 
of these sequences to those in public databases was obtained 
using Blast analysis (http://www.ncbi.nlm.nih.gov/BLAST/, last 
accessed January 24, 2014). 

Methylation-Sensitive Amplification 
Polymorphism Analysis 

Alterations in DNA methylation were investigated by MSAP 
analysis (Xiong et al. 1999; Benhattar and Clement 2004; 
Roberts et al. 2007). Aliquots of 500 ng genomic DNA of 
the parental, amphihaploid, and amphidiploid plants were di- 
gested with 1 0 U EcoRI and 20 U HpaW or 1 0 U EcoRI and 1 0 U 
Mspl (NEB) at 37 °C for 1 2 h, and the digested fragments were 
ligated to 5 pmol EcoRI and 50 pmol Hpa\\/Msp\ adaptors (sup- 
plementary table S1, Supplementary Material online), as de- 
scribed for the AFLP method. The amplicons amplified from 
this template using EcoRI and Hpa\\/Msp\ primers lacking any 
selective bases (supplementary table S1, Supplementary 
Material online) were diluted 1:30 in ddH 2 0 and used as the 
template for subsequent selective amplifications, for which 
fluorescently labeled EcoRI and unlabeled Hpa\\/Msp\ primers 
each carrying three selective bases were used (supplementary 
table S1, Supplementary Material online). The nine primer 
combinations used were EcoRI selective primer #2 plus 
Hpa\\/Msp\ selective primer #5 (abbreviated "E2 + HM5"), 
E4+HM3, E4+HM7, E4+HM8, E6 + HM6, E6+HM8, 
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E7 + HM1, E8+HM2, and E8+HM8. Electrophoresis and 
fragment detection were carried out in the same way as for 
the AFLP analysis, and each assay was repeated twice. For 
statistical analysis of the differences between the midparental 
value and the hybrids (U values), the following formula was 
used (Zhao et al. 2007): 

yi+J^ } _ s = pq(— + — ) ; 
,,_Pi - Pi 

°P1 -P2 

where n-i represents the total number of fragments of the 
midparent values, n 2 the total fragments of a hybrid, y<i the 



total number of methylated and hemimethylated fragments 
of the midparent values, y 2 the total number of methylated 
and hemimethylated fragments in the hybrid, p<\ the propor- 
tion (%) of methylated and hemimethylated fragments of the 
midparent values, and p 2 the proportion (%) of methylated 
and hemimethylated fragments in the hybrid. 

Results 

Wide Cross and Ovule Rescue 

By 12 h after pollination, pollen germination had been initi- 
ated, and the L. paludosum pollen tubes had begun to pen- 
etrate the C. morifolium stigma (fig. 2A). The hybrid zygote 
developed normally up to the multicell proembryo stage by 1 0 




Fig. 2. — The germination of L. paludosum pollen on the Chrysanthemum morifolium stigma and the development of the hybrid embryo. (A) Pollen 
germination 1 2 h postpollination. Bar: 1 00 \im. (B, Q Two globular stage proembryos with a degenerating endosperm in a 1 2-DAP ovule. (D) A degenerated 
globular proembryo in an 1 8-DAP ovule. Bar: 20 (im. (E-H) The mitotic chromosomes of (£) C. morifolium cv. "Zhongshanzigui" {In = 54), (/=) Leucanthemum 
paludosum (2n = 1 8), (G) the C. morifolium x L. paludosum amphihaploid {In = 36), and (/-/) the C. morifolium x L. paludosum amphidiploid (2n = 72). Bar: 
5 |im. (/) GISH analysis of the C. morifolium x L. paludosum amphihaploid. The chromosomes were probed with labeled C. zawadskii genomic DNA; the 27 
yellow fluorescing chromosomes were inherited from C. morifolium, and the nine unlabeled (red) ones from L. paludosum. Bar: 5|im. 
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Table 1 



Nine Annotated Novel AFLP Fragments Affected by Either Hybridization or Polyploidization 


Fragment Size (bp) 


Fragment Display Pattern 




Sequence Similarity 3 


C m 


Hy 


Am 


L.p 




189 


+ 


+ 




Ricinus communis hypothetical protein, mRNA 


122 


+ 






Acrossocheilus monticola microsatellite sequence 


124 + b 


+ 




+ b 


Helianthus annuus NBS-LRR resistance-like protein 


188 + b 




+ 


+ b 


Linum usitatissimum microsatellite sequence 


173 




+ 




Linum usitatissimum microsatellite sequence 


92 c + b 


+ 




+ b 


Linum usitatissimum microsatellite sequence 


125 




+ 




PREDICTED; Strongylocentrotus purpuratus uncharacterized, mRNA 


211 + b 


+ 




+ b 


Linum usitatissimum microsatellite sequence 


188 + b 




+ 


+ b 


Linum usitatissimum microsatellite sequence 



Note. — C. m, C morifolium; Hy, amphihaploid; Am, amphidiploid; L. p, L. paludosum; +, parental fragment present; -, parental fragment absent. Detailed information 
displayed in supplementary table S2, Supplementary Material online. 
a BlastN or BlastX threshold E values <0.05. 
b Nonparental fragments. 

c Appeared twice in different amphihaploid lines. 



DAP. However, by 1 2 DAP, only two (0.5%) out of 400 ovules 
formed globular proembryos (fig. 2B and 0- No further de- 
velopment of embryos occurred beyond this time point, after 
which embryo abortion occurred completely (fig. 2D). As a 
result, to obtain amphihaploid plants, ovules were rescued 
at 12 DAP. Of the 1,287 ovules cultured in vitro, only 105 
(8.2%) formed callus, and of these only three (2.9%) gener- 
ated buds and roots. The mature plants of amphihaploid 
and amphidiploid were morphologically homogeneous 
but were distinct from either of the parents (fig. 1). The 
plant height, plant canopy, flower diameter, and leaf size 
of the amphihaploid were smaller than those of the two pa- 
rental plants. The plant height, plant canopy, and flower di- 
ameter were improved in amphidiploid compared with the 
amphihaploid, while those were still inferior to the two paren- 
tal plants. 

Chromosome Number and GISH Analysis 

The somatic chromosome number of C. morifolium is 54 
(fig. 2£), whereas that of L paludosum is 18 (fig. 2F). As ex- 
pected, the amphihaploid and amphidiploid somatic chromo- 
some numbers were, respectively, 36 and 72 (fig. 2G and H). 
GISH analysis confirmed the hybridity of the three puta- 
tive amphihaploids; when probed with labeled C. zawadskii 
genomic DNA, 27 of the mitotic chromosomes from 
"Zhongshanzigui" recognized the labeled probes, whereas 
the other 9 from L paludosum did not (fig. 21). There was 
no evidence for any chromosome elimination. 

Genomic Changes Detected under Intergeneric 
Hybridization and Polyploidization Cycle 

All three amphihaploids appeared to be cytologically uniform 
and shared a similar plant phenotype; hence, the investigation 



into genomic perturbation was conducted on these three 
plants. The nine AFLP primer combinations amplified 375 frag- 
ments from the C. morifolium and 354 from the L. paludosum 
template. The AFLP profiles of amphihaploid lines 1/2/3 con- 
sisted of 388/396/389 fragments, respectively, of which 4/3/3 
were not present in either parent's profile (table 1, supple- 
mentary table S2, Supplementary Material online), whereas 
118/115/115 (30.4%/29.0%/29.6%) were present in both 
parental profiles (fig. 3AM 180/188/182 (46.4%/47.5%/ 
46.8%) were inherited from C. morifolium (fig. 3A2), 
and 86/90/89 (22.2%/22.7%/22.9%) from L paludosum 
(fig. 3/43). The amphidiploid's AFLP profile comprised 
389/396/390 fragments, 381/391/384 of which were in- 
herited from one of the parental lines. Of the 4/3/3 nonpa- 
rental fragments detected in the amphihaploid, 3/3/3 were 
also present in the amphidiploid's profile, but only one was 
absent in the amphidiploid line 1 . The 3/2/2 parental frag- 
ments present in the amphihaploid's profile were absent 
from the amphidiploid's profile, whereas a further 3/1/2 pa- 
ternal parent fragments were present in the amphidiploid's 
profile but not in the amphihaploid's. A further 2/1/1 frag- 
ments were unique to the amphidiploid's profile (fig. 3D). Of 
the 27 nonparental fragments, two fragments yielded poor 
sequence quality and 25 fragments were successfully se- 
quenced. Fifteen out of 25 fragments showed no homology 
to any GenBank entries, two matched hypothetical or unchar- 
acterized genes, and the remaining eight sequences hit micro- 
satellite sequence and a sequence encoding an NBS-LRR 
resistance-like protein (table 1, supplementary table S2, 
Supplementary Material online). Remarkably, three fragments 
were present twice in different amphihaploid lines, that is, one 
fragment was present in lines 1 and 3, and the remaining two 
were present in lines 2 and 3 (supplementary table S2, 
Supplementary Material online). 
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AFLP SRAP SCoT 




Fig. 3. — Representative DNA profiling of the Chrysanthemum morifolium x Leucanthemum paludosum amphihaploid, amphidiploid, and its parents. 
(A) AFLP. (B) SRAP, (0 ScoT: 1, fragments present in both parental profiles; 2, fragments inherited from C. morifolium; and 3, fragments inherited from L 
paludosum. (D) Variation in AFLP profiles, from top to bottom: most of the parental fragments were inherited by both the amphihaploid and the amphi- 
diploid; four nonparental fragments present in the amphihaploid were also present in the amphidiploid; a nonparental fragment detected in the amphi- 
haploid but not transmitted to the amphidiploid; three parental fragments present in the amphihaploids but not in the amphidiploid; three parental 
fragments absent in the amphihaploid but present in the amphidiploid; and two nonparental fragments present in the amphidiploid but not in the 
amphihaploid. (E) Variation in SRAP profiles (line 1 ), from top to bottom: two nonparental fragments present in both the amphihaploid and the amphidiploid; 
two parental fragments not present in either the amphihaploid or the amphidiploid; and one nonparental fragment present in the amphidiploid but not in 
the amphihaploid. 



Table 2 



SRAP Fragments Affected by Either Hybridization or Polyploidization in Amphihaploid Line 1 


Fragment Size (bp) 


Fragment Display Pattern 




Sequence Similarity 3 


C m 


Hy 


Am 


Lp 




216 


+ 


+ 




PREDICTED: Ciona intestinalis similar to HyTSRI protein mRNA 


125 


+ 


+ 




None 


138 + b 


+ 




+ b 


None 


131 + b 


+ 






Cynara cardunculus cardosin B gene for aspartic proteinase, exons 1-14 


160 




+ 




Brassica rapa subsp. pekinensis clone KBrB084K02, complete sequence 



Note. — C. m, C morifolium; Hy, amphihaploid; Am, amphidiploid; L p, L paludosum; +, parental fragment present; -, parental fragment absent. Detailed information 
displayed in supplementary table S2, Supplementary Material online. 
a BlastN or BlastX threshold E values <0.05. 
b Nonparental fragments 



The 24 SRAP primer combinations amplified 318/322/326 
fragments from the three amphihaploid templates, of which 
101/103/106 (31.8%/32.0%/32.5%) were present in both 
parents' profiles (fig. 3£1), 120/121/117 (37.7%/37.6%/ 
35.9%) were inherited from C. morifolium (fig. 3B2), 
and 95/95/101 (29.9%/29.1 %/31 .0%) from L paludosum 
(fig. 3B3). The 2/3/2 (0.6%/0.9%/0.6%) novel fragments 
were detected (table 2, supplementary table S2, 
Supplementary Material online), all of which were transmitted 



to the amphidiploid. The amphiploid's profiles also showed 
the loss of 2/1/1 parental fragments and the gain of 1/1/1 
novel fragment. However, only a few fragments can find se- 
quence similarity with accessions deposited in the GenBank, 
and no repeat sequences between amplified fragments have 
been found (supplementary table S2, Supplementary Material 
online). For example, in amphihaploid line 1, only three 
nonparental fragments showed homology to Ciona intestina- 
lis HyTSRI protein mRNA, Cynara cardunculus cardosin B gene 
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for aspartic proteinase, and a clone of Brassica rapa subsp. 
pekinensis sequence (fig. 3E and table 2). The ScoT profiles 
(12 primers) comprised 200/210/202 fragments (the amphi- 
haploid) and 199/207/202 (the amphidiploid). These frag- 
ments were classifiable as 1/2/2 (0.5%/1.0%/1.0%) 
nonparental, 65/70/67 (32.5%/33.3%/33.2%) inherited 
from either parent (fig. 3CI), 77/80/75 (38.5%/38.1 %/ 
37.1%) from C. morifolium (fig. 3C2), and 57/58/58 
(28.5%/27.6%/28.7%) from L paludosum (fig. 3C3); the 
difference in fragment number between the amphihaploid 
and amphidiploid profiles showed a loss of 2/4/2 and a gain 
of 1/1/2 fragments in amphidiploid. 

Epigenetic Changes Induced by Allopolyploidization 

Changes in methylation status were investigated using MSAP 
profiling, based on nine combinations of EcoRI + Hpa\\/Msp\ 
primer combinations (fig. 4). Four types of fragment were 

C. morifolium amphihaploid amphidiploid L paludosum 




Non-methylated 
Fullymethylated 
Hemimethylated 
Hypermethylated* 



Fig. 4. — Representative variation in MSAP profiles. " Variation of 
DNA methylation between amphihaploid and amphidiploid; "+," frag- 
ments obtained after digestion with EcoR\ or /-/pall. fragments not 
digested by EcoRI or /-/pall. Note: hypermethylation (type IV) is detectable 
from the comparison between the amphihaploid and amphidiploid pro- 
files. For example, fragments absent from amphihaploid (both H and M 
lanes) but present in amphidiploid (either H or M lane), and vice versa. 



revealed: type I fragments are nonmethylated and were pre- 
sent in both the H and M digests; type II are fully methylated 
and only appeared in the M digests; type III are hemimethyl- 
ated and appeared in the H digests; finally, type IV were frag- 
ments absent from both H and M lanes in amphihaploid but 
present in either H or M lane of amphidiploid, inferring de- 
creased methylation, or those absent from both H and M 
lanes in amphidiploid but present in either H or M lane of 
amphihaploid, indicating increased methylation. 

The C. morifolium profile included 267 type I, 169 type II, 
and 1 14 type III fragments, whereas the corresponding num- 
bers in L. paludosum were 264, 1 54, and 1 1 6. The proportion 
of methylated fragments was 51.5% in C. morifolium and 
50.6% in L. paludosum. In the amphihaploid lines, the propor- 
tion of cytosine methylated fragments was 45.1-45.2%, a 
level that is inconsistent with the midparent value (table 3). 
The amphihaploids included a significantly lower proportion of 
fully methylated fragments (23.4-23.8%) than the midparent 
level (the resulting U values were 2.02-2.23, with 
L/ 0 .o5= 1.96); however, with respect to the proportion 
of hemimethylated fragments present, the proportion of 
amphihaploid templates (20.8-21.3%, with a U value was 
0.01-0.13) equaled the midparent value. The mean number 
of fragments that showed changes in methylation state be- 
tween the amphihaploid to amphidiploid profiles was 161.3 
on average, of which 57.3 showed a decrease in methylation, 
whereas the remainder exhibited an increase in methylation 
(table 3). It showed that 21.3 of these sites involved a shift 
from type III to type I and 1 4.3 from type IV to type II (decreased 
methylation), whereas 39.0 shifted from type I to type III and 
23.0 from type I to type II (increased methylation) (table 4). 

Discussion 

Intergeneric Hybridization between C morifolium and 
L paludosum 

Embryo rescue is a frequently used technique in wide cross, 
which has allowed for the successful production of many 
interspecific and intergeneric hybrids (Van Tuyl et al. 1991; 
Sharma 1995; Tang et al. 2009). Although many 
C. morifolium x L paludosum hybrid zygotes were able to 



Table 3 



Levels of Cytosine Methylation Detected in the Amphihaploid and Parental Plants 


Plant Lines 


Total Sites 


Nonmethylated 




Methylated 








Type I 


Total (II + III) 


Type II 


Type III 


C m 


550 


267 (48.5%) 


283 (51.5%) 


169 (30.8%) 


114 (20.7%) 


Lp 


534 


264 (49.4%) 


270 (50.6%) 


154 (28.8%) 


116 (21.7%) 


Midparent value 


100% 


49.0% 


51.0% 


29.8% 


21.2% 


Hy1 


554 


304 (54.9%) 


250 (45.1%) 


132 (23.8%) a 


118 (21.3%) 


Hy2 


556 


306 (55.0%) 


251 (45.0%) 


136 (24.5%) a 


114 (20.5%) 


Hy3 


575 


311(54.1%) 


264 (45.9%) 


141 (24.5%) a 


123 (21.4%) 



a Values significantly lower than the corresponding midparent value at the U 0 _ 05 threshold. 
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Table 4 

MSAP Fragments Affected by Either Hybridization or Polyploidization 



Fragment Type 




Fragment Display Pattern in MSAP Gel 




Number of Sites 


Status 


Hy 


Am 


Hy H Lane 


Hy M Lane 


AM H Lane 


AM M Lane 


Type IV 


Type 1 






+ 


+ 


5.0 ±0.9 




Type 1 


Type IV 


+ 


+ 






14.7 ±0.5 


t 


Type IV 


Type III 






+ 




1 1 .0 ± 1 .4 


I 


Type III 


Type IV 


+ 








20.7 ±1.0 


t 


Type IV 


Type II 








+ 


14.3±1.2 


I 


Type II 


Type IV 




+ 






7.0 ±0.9 


t 


Type II 


Type 1 




+ 


+ 


+ 


5.3 ±0.3 


I 


Type 1 


Type II 


+ 


+ 




+ 


23.0 ±1.9 


t 


Type III 


Type 1 


+ 




+ 


+ 


21.3±1.7 




Type 1 


Type III 


+ 


+ 


+ 




39.0 ±2.8 


t 



Note. — |, decreased methylation; f, increased methylation; Hy, amphihaploid; Am, amphidiploid. 



develop the multicellular proembryo, further development of 
embryo did not occur in vivo. Therefore, we resorted to ovary 
culture, which contributed to successful hybrid zygotes rescue 
before hybrid embryo abortion occurred. However, only a small 
number of true hybrids could be regenerated to form viable 
amphihaploid seedlings. The amphihaploid and amphidiploid 
are able to flower normally, but the meiotic pairing behavior of 
the chromosomes is disordered in the amphihaploid. 

Genetic Relationship between Chrysanthemum and 
Leucanthemum 

The meiotic pairing behavior of the chromosomes in a wide 
hybrid is arguably the most reliable means to deduce the tax- 
onomic relationship between the parental species; however, 
the taxonomic relationships remain as yet not fully resolved 
within Chrysanthemum sensu lato (Teixeira da Silva 2003). In 
the absence of such data, GISH analysis can provide a strong 
basis for assessing the level of distinctness between pairs of 
genomes, and it has been widely used to characterize hybrids 
generated from species of the genera Chrysanthemum, 
Ajania, Tanacetum, Brachanthemum, Elachanthemum, 
Leucanthemella, and Nipponanthemum (El-Twab et al. 
1999; El-Twab and Kondo 2001; Zhao et al. 2008). 
Genomes that are closely related to one another require the 
incorporation of blocking (unlabeled) DNA to achieve success- 
ful discrimination (El-Twab and Kondo 2004; Deng et al. 
2010). Here, GISH, in which the probe DNA was obtained 
from the diploid C. zawadskii, was able to distinguish the 27 
C. morifolium chromosomes from the 9 L paludosum ones in 
mitotic spreads of the amphihaploid plant without any block- 
ing DNA. Chrysanthemum zawadskii has been proposed as a 
potential ancestor of C. morifolium (Yang et al. 2006); here, 
our data again suggested a relatively close relationship be- 
tween the two species within Chrysanthemum, whereas the 
relationship between Chrysanthemum and Leucanthemum 
spp. in the tribe Anthemideae is distant, which is consistent 
with the observation of hybrid embryo abortion. 



Genome Perturbation in the Synthetic Amphihaploids 
and Amphidiploids 

Plants tolerate hybridization and polyploidization much more 
easily than do animals. Newly synthesized plant polyploids 
have been observed to undergo a range of chromosomal rear- 
rangements, the activation of dormant transposable elements, 
the elimination of DNA sequences, and epigenetic silencing 
(Barton 2001; Chen 2007; Xiong et al. 2011). Commonly, 
wide hybrids express certain characteristics (at either the mor- 
phological or molecular level) clearly inherited from one of the 
parents, some from the other parent, and a number of de 
novo ones that are apparently not inherited from either 
parent (Soltis and Soltis 1999; Richards 2003; Bell et al. 
2013). Synthetic allopolyploids between wheat and its near 
relatives suffer a measurable degree of sequence elimination 
(Shaked et al. 2001), as also do those generated in Cucumis 
spp. (Chen et al. 2007), Brassica spp. (Song et al. 1995), and 
Tragopogon spp. (Tate et al. 2006). Here, AFLP profiling 
showed that DNA diversification had probably occurred in 
the C. morifolium x L paludosum amphidiploid, as some frag- 
ments were lost, whereas others were gained in amphidiploid. 
Several novel fragments were also detected in the amphiha- 
ploid, showing that genome changes are also induced by the 
intergeneric hybridization event itself. New fragments can also 
form via point mutation if the mutations remove (or create) a 
restriction site. Genomic shock can activate quiescent transpo- 
sons/retrotransposon and some mutation (Liu et al. 1998; 
Shaked et al. 2001; Kashkush et al. 2002). We did not find 
any new fragments corresponding to retrotransposon or DNA 
transposons, however, some fragments contained microsatel- 
lite sequence, suggesting whole-genome doubling may be 
accompanied by the modification of microsatellite sequence. 

Although marker-based analysis has become the method 
for many studies on genetics, evolution, and ecology, 
genotyping error rates are not explicitly calculated in most 
empirical studies (Pompanon et al. 2005). These errors cause 
confusion in the signal-to-noise ratio and hence a loss of 
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resolving power; fortunately, they might not bias all the results 
of the analysis (Bonin et al. 2004). The number of amplifica- 
tion fragments in each primer combination is decided by 
the template quality, primer sequence, PCR reaction condi- 
tions, competition between potential amplicons, and base 
mismatching between primer and template (Smith and 
Williams 1994; Bussell et al. 2005). In this study, we attempt 
to minimize error rates by employing optimized DNA extrac- 
tion, reaction mix, PCR reaction conditions, and primer screen- 
ing. Theoretically, the rate of genetic diversity for genotypes 
was nearly the same using AFLP, SRAP, or SCoT markers 
(Schuler et al. 1996; Arabidopsis Genome Initiative 2000; Li 
and Quiros 2001 ; Meudt and Clarke 2007). In this study, how- 
ever, the proportion of each parent's contribution to the 
amphihaploid profiles differed when AFLP and SRAP or ScoT 
markers were used, and only 46.4% of the amphihaploid's 
AFLP fragments were inherited from C. morifolium, whereas 
22.2% were inherited from L paludosum; in contrast, the 
relative contributions with respect to the SRAP fragments 
were 37.7% and 29.9%, and for the SCoT fragments were 
38.5% and 28.5%. The difference between these contribu- 
tions implies a difference in the nature of the fragments de- 
tected by the three profiling methods; a plausible explanation 
could relate to the random nature of AFLP fragments, which 
reflect genome-wide variation in DNA sequence (largely con- 
sisting of noncoding DNA), whereas both SRAP and SCoT 
fragments are derived mostly from coding sequence. The 
reason may be that some "redundant" sequence possibly 
from the noncoding region of L paludosum was eliminated 
in the amphihaploid and the change may be nonrandom 
during intergeneric hybridization. Sequence elimination is 
one of the earliest responses of the genome to wide hybrid- 
ization or allopolyploidy (Shaked et al. 2001). The identity of 
the sequences eliminated has been shown in wheat (Shaked 
et al. 2001; Kashkush et al. 2002), Tragopogon spp. (Tate 
et al. 2006; Koh et al. 2010), Cucumis spp. (Chen et al. 
2007), and Brassica spp. (Song et al. 1995). Our present 
data support previous observation that sequence elimination 
after wide hybrid is reproducible (McClintock 1984). 

Cytosine methylation is the most common covalent modi- 
fication of DNA in eukaryotes, and it provides a ubiquitous 
epigenetic control (Zilberman et al. 2006). In this study, 
Chrysanthemum and Leucanthemum methylated sites are 
up to 51.6% and 50.6%, respectively (table 3). A number 
of epigenetic phenomena have been associated with wide 
hybridization and whole-genome doubling (Lee and Chen 
2001; Wang et al. 2006). For example, >13% of the site 
alterations in cytosine methylation occurred either in the 
hybrid or in the allopolyploid in wheat (Shaked et al. 2001). 
After hybridization, the global level of DNA methylation, as 
measured by MSAP, was somewhat lower in the amphiha- 
ploid than in the parental genome, which suggests that one 
effect of the wide hybridization may be to reduce global 
DNA methylation. Similar epigenetic phenomena have been 



observed in other newly synthesized triticale (wheat x rye 
hybrids), which tended to be rather unstable in the first few 
generations, probably owing to a reduced methylation in the 
hybrids, but eventually became stable (Ma and Gustafson 
2008). It should be noted that a drop in DNA methylation 
could also be induced by the callus phase that the primary 
hybrids go through (Miguel and Marum 2011). However, it 
has been reported that in spite of the different kind of meth- 
ylation changes in individual regenerants of Codonopsis lan- 
ceolata, overall tissue culture did not cause significant global 
alteration in the DNA methylation (Guo et al. 2007). In oil 
palm, the highest number of polymorphic bands had occurred 
during tissue culture, and most frequently, the HpaW recogni- 
tion was only 0.3% (Matthes et al. 2001). Further investiga- 
tion using whole-genome bisulfite sequencing revealed that 
these changes resulting from tissue culture accounted for a 
relatively small proportion of the rice genome (Stroud et al. 
2013). We suggested that wide hybridization may be the 
major driver of global DNA methylation reduction. 

DNA methylation is generally recognized to suppress gene 
expression as regulatory factors, and changes in cytosine 
methylation resulting from wide hybridization and/or polyploi- 
dization can be expected to have a major effect on gene ex- 
pression in polyploids. It has long been recognized that genes 
can be silenced by methylation and turned on by demethyla- 
tion (Martienssen and Colot 2001). Here, the global pattern of 
methylation in the amphihaploid was noticeably different than 
those of the two parental lines (table 3). Although it was 
reduced by the wide hybridization event, it was restored 
somewhat by whole-genome doubling, and the fragments 
of increased DNA methylation in the amphihaploid are 1 .8 
(104/57.3) times those of decreased methylation. Similarly, 
in the allotetraploid of Cucumis x hytivus, 68.2% (15 from 
22 sites) of methylation loci changes showed an increase in 
cytosine methylation, whereas only 31 .8% displayed demeth- 
ylation (Belzile et al. 2008). Taken together, it suggests that 
DNA methylation partly functions epigenetically and dynami- 
cally to control and compromise with unbalanced gene 
expressions caused under certain circumstances. Although 
the underlying mechanisms for these phenomena are not 
yet determined, it is possible that the cytosine methylation 
machinery responds as a "genome defense" system by unit- 
ing divergent genomes into one nucleus (Yoder et al. 1997; 
Osborn et al. 2003). 

Conclusion 

The taxonomic relationship between Chrysanthemum and 
Leucanthemum is relatively distant. Combining such diverse 
genomes into a single nucleus can be expected to cause 
genomic shock in the sense of McClintock: "a highly pro- 
grammed sequence of events within the cell that serves to 
cushion the effects of the shock." The induced changes to 
the genome must occur very early in the development of the 
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hybrid zygote and continue when the amphihaploid is sub- 
jected to whole-genome doubling. A further manifestation of 
the shock is the alteration in the pattern of DNA methylation. 
These rapid and massive changes in the genome of hybrids 
may serve as a tolerance of genome shock and might contrib- 
ute to rapid genetic diploidization of the newly formed allo- 
polyploid. An important question is whether the genomic and 
epigenomic events observed here are common in wide hybrids 
within the Asteraceae family. Understanding these genomic 
changes and epigenetic modifications is a worthy goal to 
increase understanding of the evolutionary process of the 
Asteraceae family, which will require the production and anal- 
ysis of other wide hybridization combinations between mem- 
bers of the Asteraceae family. 

Supplementary Material 

Supplementary tables S1 and S2 are available at Genome 
Biology and Evolution online (http://www.gbe.oxfordjournals. 
org/). 
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